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ABSTRACT 

It has been suggested that radiative transfer effects may explain the unusually high equivalent widths 
(EWs) of the Lya line, observed occasionally from starburst galaxies, especially at high redshifts. If 
the dust is locked up inside high-density clouds dispersed in an empty intercloud medium, the Lya 
photons could scatter off of the surfaces of the clouds, effectively having their journey confined to the 
dustlcss medium. The continuum radiation, on the other hand, does not scatter, and would thus be 
subject to absorption inside the clouds. 

This scenario is routinely invoked when Lya EWs higher than what is expected theoretically are 
observed, although the ideal conditions under which the results are derived usually are not considered. 
Here we systematically examine the relevant physical parameters in this idealized framework, testing 
whether any astrophysically realistic scenarios may lead to such an effect. It is found that although 
dumpiness indeed facilitates the escape of Lya, it is highly unlikely that any real interstellar media 
should result in a preferential escape of Lya over continuum radiation. Other possible causes are 
discussed, and it is concluded that the observed high EWs are more likely to be caused by cooling 
radiation from cold accretion and/or anisotropic escape of the Lya radiation. 
Subject headings: radiative transfer — scattering — galaxies: ISM 



1. INTRODUCTION 

Our understanding of the high-redshift Universe has 
expanded tremendously during the last decade, espe- 
cially due a single specific emission line, the Lya line. 
Tracing in particular galaxies in the process of forming, 
the shape, intensity, and spatial distribution of this line 
carry a wealth of information. However, while our theo- 
retical understanding of the physical processes that gov- 
ern the radiative transfer (RT) of Lya has also advanced 
significantly, it is evident that the complexity of density 
fields, gas kinematics, dust distribution, line formation, 
impact of the intergalactic medium (IGM), etc. render 
individually observed lines quite difficult to interpret in 
detail. 

As Lya emitt ing galaxies (LAE s) generally are 
relatively small (iGawiser et all 120061: lLai et all 120071: 
INilsson et al.1 12007 ) of low (specific) star formation 
rate (IFvnbo et al.1 120031: INilsson et all 120101 ) . our lack 
of knowledge is in part due to the difficulties in 
observing faint galaxies. This obstacle may to 
some extend be overcome using g ravitational!/ l ensed 
LAEs (e.g. iFosburv et all [20031: lOuider et all [200l 
iChristensen et all l2012al lbf). which may boost the ob- 
served flux by more than an order of magnitude. Larger 
surveys of the statistical properti es of LAEs, such 
the ir luminosity fu n ction (LF. lOuchi et aTl 
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our ignorance of RT processes easily may result in sys- 
tematic errors. A common technique by which to as- 
sess the stellar population of individual galaxies at high 
redshifts is to measu re the equivalent width (EW) W 
of the Lya line (e.g. IGronwall et all 120071: lOuchi et all 
120081: iStark et al.il2010l) . This quantity is defined as the 
ratio of the integrated line flux to continuum flux den- 
sity, W = JdX [F(X) — F a ]/F , where the exact limits 
of the integral are not important, as long as the full 
line is included. For sloped continua, the wavelength 
dependency of the continuum must also be taken into 
account. For high EWs, W can be approximated by the 
relative escape fraction of Lya and continuum photons. 
The EW depends on galactic parameters such as the ini- 
tial mass func tion (IMF) and the gas metallicity (e.g. 
ISchaererl l2002). and can consequently be used as a probe 
of these quantities. 

In young galaxies, Lya is produced mainly from re- 
combinations following the ionization of hydrogen sur- 
rounding O and B stars. Since these stars are short- 
lived, a few Myr after an initial starburst the EW 
will generally decrease significantly . Stellar population 
syntheses (e.g. ICharlot fc FaTll fl993l: IValls-Gabaua1ll99l 
ISchaererl 120031 ) predict that the EW should initially 
reach ~ 240 A, eventually declining and settling on 
roughly 80 A, where the exact values depend on the 
assumed IMF and metallicity of the population. Cu- 
riously, several observational studies have reported on 
the detection of much higher EWs (e.g. IHu k, McMaho"nl 
19961: iKuTritzki et all [2000t IMalhotra fc Rhoadsl 120021: 
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While some of these extreme EW galaxies can be at- 
tributed to AGN activity, in many cases this explanation 
is explicitly excluded. Moreover, such high EWs are not 
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readily explained by simply assuming a more top-heavy 
IMF. Indeed, the observations pose a serious challenge 
to our understanding of both galaxy formation, stellar 
evolution, and radiative transfer. 

In this work, we will distinguish between three related 
notations. The equivalent width W m t of the intrinsically 
emitted spectrum is given by the stellar population; for 
instance, a more top-heavy IMF emits a harder spec- 
trum, with relatively more ionizing and hence Lya pho- 
tons. When escaping the galaxy, the spectrum has an 
equivalent width Wemi as we shall see below, the ISM 
may possess the ability of transferring an uneven frac- 
tion of Lya vs. continuum photons, making W em larger 
or smaller than W- m t- Finally, as the total number of 
photons in a line is conserved when traveling through 
the Universe, whereas the number of continuum photons 
per wavelength bin is reduced as they are cosmologically 
rcdshiftcd, the observed equivalent width Wobs from a 
source at redshift z is given by W bs = (1 + z)W cm . 

1.1. The Neufeld scenario 

As Lya photons scatter on neutral hydrogen, their 
path length before escaping a galaxy will always be longer 
than that of continuum radiation. Thus, the immediate 
corollary is that Lya will be more susceptible to dust ab- 
sorption thaiithe continuum, implying that W om < W- ln t. 
However . iNeufeldl ([ 19911 . hereafter N91) investigated ana- 
lytically the resonant scattering of Lya photons through 
the ISM and found that under special circumstances, the 
Lya photon may actually suffer less attenuation than 
radiation which is no t resonantly scattered, e.g. contin- 
uum radiation. IN91I considered the escape of radiation 
from the center of a plane-parallel, two-phase structure, 
in which sperical clouds of homogeneously mixed neu- 
tral hydrogen and dust lie embedded within a virtually 
empty intercloud medium (ICM). The number density of 
clouds arc assumed to be small enough that they do not 
touch, yet sufficiently numerous that they cover most of 
the sky. 

Furthermore, the source of both Lya and continuum 
photons is assumed to be pointlike and situated in the 
cent er of the slab, in the ICM. Under these conditions 
IN91I showed that Lya photons, upon entering a cloud, 
scatters only a few times before returning to the ICM, 
thereby being exposed arbitrarily little to absorption by 
dust. Consequently, the journey of the Lya photons will 
be confined primarily to the dustless ICM, preserving the 
total Lya luminosity. In contrast, the continuum radia- 
tion which penetrates the clouds rather than scattering 
off of their surfaces will be subject to the full attenua- 
tion, given by the optical depth of dust from the center 
and out. 

Thus, since the continuum is reduced while the line is 
more or less preserved, the EW of the escaping radiation 
can be "boosted" to arbitrarily high values. We define 
such a boost as 

b^^. (1) 

Of course, the boost does not imply an increased Lya 
luminosity, but rather a reduced continuum luminosity. 

1.2. Numerical and observational support 



The m ere poten tial of an EW boost w as the conclusion 
of Hp's studies. lHansen fc Ohl (pOOl . hereafter HO06) 
studied the scenario in more detail, both analytically and 
numerically. In addition to confirming the boost, they 
also looked into the effect of galactic outflows of clumpy 
gas on the line profile, exploring the RT in various ge- 
ometries; 

While [HO06I did investigate various extensions of the 
Neufeld scena rio, one of their main conclusion was the 
same as IN91I : "If most of the dust resides in a neu- 
tral phase which is optically thick to Lya, the Lya 
EW can b e strongly enhanced". The results of IN91I 
and IHO06I seem a natural explanation of the unusually 
high observed EWs, and lately such observations has 
routinely been interpreted within the framework of this 
scenario, inferring the pres ence, or non-pr e sence, of a 
clump y, dusty ISM Thu s, IChapman et all ([20051 ) and 
lately [Bridge et all (|2012[ ). interpret the fact that Lya 
radiation is visible from ~ 50% of a sample of sub- 
mm galaxies as suggesting the presen ce of very patchy 
and in homogeneous dust distribution. iFinkelstein et aLl 
(|2007[) regards a clumpy ISM as a possible explanation 
of their high Lya EWs, and in IFinkelstein et al.1 ([20081 ) 
propose a parameter q to characterize dumpiness, such 
that Lya is attenuated by a factor e~ qT , where r is 
the dust optical depth for non-resonant photons close 
to the Lya line. For a given Lya EW observation 
they then calculate the dumpiness of the ISM, positing 
that for q < 1 the EW is boosted du e to dust dumpi- 
ness. This method is furthe r used in IFinkelstein et all 
(l2009aB.[20TIal) . IYuma eTall (l20Tol).IBlanc et all (12011D . 
iNakaiima et all (|2012l) and lHashimoto et al.l (|2012f ) to 
infer various degrees of dumpiness in observed LAEs. A 
equi valent approach w as followed in iNiino et al.l (12 009) 
and iKobayashi et al.1 ([20101 ) . Similarly. iDavkl et al.l 
([20081120091 12010L 120111 ) lias invoked a clumpy ISM sce- 
nario in order to match their simulated LAE LFs to ob- 
served ones. 

However, the idealized conditions under which a boost 
may be achieved is rarely taken into account. Accord- 
ingly, we believe that a revision of the scenario is timely, 
and in the present work we aim to investigate the valid- 
ity of the various premises of the model. We do this by 
calculating numerically the RT of Lya and continuum 
in model galaxies constructed to cover a broad range 
of plausible as well as i mpla usible systems. In addi- 
tion to the work of lN9ll and IHO06I . the topic of reso- 
nant scattering in a multiphase medium has received at- 
tention from iRichling et all (120031). iSurlan et all (I2012D . 
iDiikstra fc Kramerl (I2012D. an d lDuval et al.l (|2012l ). Of 
these. onlv lDuval et al.l (|2012l) concern themselves with 
the relative escape fraction of Lya and continuum, and 
found tentative evidence for the impracticability of en- 
hancing the Lya EW, by investigating the Lya and con- 
tinuum RT in an expanding shell of clumps with different 
values of dust optical depth, expansion velocity, Hi den- 
sity, and line width. They found that indeed a virtually 
empty ICM was needed, along with only modest outflow 
velocities and high E(B — V), for an ISM geometry to 
result in a boost. 

After describing the principles of the applied numeri- 
cal code in Sec. [21 we systematically vary all parameters 
relevant for decribing a (model) galaxy in Sec. [3j one at 
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a time. This gives us a feeling for the impact a given pa- 
rameter has on the boost. However, as the change of one 
parameter may very well be either enhanced or counter- 
acted by the change of another, after discussing in Sec. 0] 
observational and theoretical constraints on the actual 
values of the parameters that are likely to be met in 
real astrophysical situations, we subsequently, in Sec. O 
undertake a large sample of simulations with random pa- 
rameters covering these values, constructing a likelihood 
map of achievable boosts. The results of these calcu- 
lations are discussed in Sec. [6l along with a discussion 
on other possible scenarios that could lead to high Lya 
EWs. 

2. RADIATIVE TRANSFER SIMULATIONS 



IN91I regarded each dusty cloud as a scattering parti- 
cle itself, with a c ertain probability of absorbing a pho- 
ton. Whcn lHO06l developed their numerical model, they 
followed the same approach. They used a Monte Carlo 
code, where the path of individual photons were traced as 
they traveled through the ICM, but whenever a cloud was 
encountered, the photon would simply be absorbed or 
scattered at once in some direction away from the cloud. 
The appropriate probability density functions (PDFs) 
were calculated on beforehand from fits to a series of 
RT calculations of photons incident on a cloud surface. 
This rendered feasible simulations that were otherwise 
impractical. In the almost seven years that have passed 
since, computational power has increased to a point that 
"brute force" simulations can in fact be performed rather 
effortlessly, obliviating the need for several approxima- 
tions. 

In the present work, the RT calculations are conducted 
using the numerical code MoCaLaTA (jLaursen et al.l 
l2009alfbT ). In the following, the basics of the code are 
outlined: 

The galaxy is constructed on a grid of cells, each of 
which holds a value of the physical parameters impor- 
tant to the RT; the neutral hydrogen density tihi, the 
dust density rid, the gas temperature T, and the three- 
dimensional velocity vt, u ik of the gas elements. The orig- 
inal Neufeld scenario concerned itself with radiation es- 
caping from a slab of gas. For numerical reasons, higher 
resolution can be achieved considering r adiation escaping 
a sphere of multiphase g did also IHO061 Qualita- 
tively, and even quantitatively except for factors of order 
unity, the results are equivalent. Thus, a galaxy is mod- 
eled as a number N c \ of spherical, non-overlapping clouds 
with radius r c i, dispersed randomly within a sphere of ra- 
dius r g ai. A cell may be either an ICM cell or a cloud cell. 
In order to make the clouds spherical, cells on the border 
between a cloud and the ICM are adaptively refined, such 
that a given cell is split into eight subcclls, recursively 
until a satisfactory resolution is achieved. The structure 
is depicted in Fig. [TJ where also a surface brightness map 
of a simulation is shown. We typically consider ~ 10 3 " 5 
clouds, with each cloud typically consisting of 10 2 " 3 cells. 
Note however that the actual shape of a cloud does not 
affect the results sign ificantl y; the most important quan- 
tities, as identified bv lHO06l . are the cloud albedo, i.e. the 
probability that a photon incident on a cloud is reflected 
rather than absorbed after a number of scatterings, and 
the average number Nq of clouds with which a Lya pho- 
ton interacts before escaping the galaxy, in the absence 




Fig. 1. — Left: Graphical representation of the design of the 
model galaxy with ~ 200 clouds. The AMR structure with six 
levels of refinement is illustrated in the inset. Note that for pre- 
sentation purposes only, a two-dimensional analogy is shown; in 
reality all simulations are carried out in three dimensions. Bight: 
Simulated surface brightness map of a model with ~ 500 rather 
large clouds of radius r c \ = 350 pc to emphasize the structure. 
The fiducial model used in Sec. \3\ has 6500 clouds of r c \ = 100 pc, 
with a maximum refinement level of ~ 10. 

of absorption. This number in turn is a function of the 
covering factor f c , which is the average number of clouds 
intercepted by a sightline from the center and out. In- 
dividual Lya and continuum photons are then emitted 
and traced as they scatter in real and frequency space 
throughout the inhomogeneous ISM. 

MoCaLaTA has previously been applied to galax- 
ies extracted from cosmological simulations. Although 
it was tes ted thoroughly agains t various analytical so- 
lutions in lLaursen et al.l (|2009al fbT) , some modifications 
and extensions had to be made in order to use it for 
the idealized simulations in the present study. These are 
described and tested in App. lAl 

For a given set of parameters N c \, r c \, etc., the observed 
boost will depend on the direction which it is observed, 
as well as on the actual random realization of cloud po- 
sitions. In addition to sampling the average spectrum, 
i.e. the spectrum of photons escaping in all directions, 
MoCaLaTA calculates the spectra escaping exactly in 
the six different directions along the Cartesian axes. The 
more clouds a galaxy comprises, the smaller the spread 
will be in different directions, and for different realiza- 
tions. Further, the variation between different directions 
is much larger than between different realizations. For in- 
stance, for a covering factor of f c = 2, only ~ 30% of the 
sightlines will actually intercept two clouds, while ~ 10% 
will intercept no clouds at all (i.e. the covering fraction 
is ~ 0.9 for f c ~ 2), and a non-negligible fraction of the 
sightlines will intercept > 5 clouds. These three cases 
will, respectively, result in a boost close to the in aver- 
age, no boost at all, and an essentially infinite boost (for 
a central point source). Thus, the fractional standard 
deviation o~b/b from different viewing angles is typically 
of order 50%, while for the 4tt average of many realiza- 
tions of the same model, cr^/b is but a few percent. In 
the following analysis, the presented values of b represent 
the Air average for a single realization of a given model. 

3. INVESTIGATING THE CRITERIA FOR A BOOST 

In order to systematize which physical conditions are 
necessary for boosting the EW, we first settle on a fidu- 
cial model. In the subsequent sections, several of the 
implied parameters are then relaxed or varied, individu- 
ally or in conjunction. When nothing else is stated, the 
remaining parameters correspond to the fiducial model. 
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Initially, we will not concern ourselves with the realism 
of the parameters, instead deferring this discussion to 
Sec. |U although we note that several of the chosen val- 
ues, as well as various derived parameters such as color 
excess and total galaxy mass, roughly resemble observed 
typical LAE values. 

The fiducial model galaxy is a sphere of radius r ga i = 5 
kpc consisting of N c \ = 6500 of clouds of equal radii 
r c i = 0.1 kpc. The density of neutral hydrogen in the 
clouds is tt-hi = 1 cm~ 3 , implying a column density 
-Whi.cI ~ 3 x 10 20 cm" 2 as measured from the center to 
the surface of a cloud. The clouds are dispersed at ran- 
dom in an empty ICM. The temperatures of the clouds 
and the ICM are T cl = 10 4 and Ti CM = 10 6 K, re- 
spectively. The physical significance of the clouds and 
the ICM are the interstellar phases conventionally called 
the warm neutral medium (WN M) and the hot ion ized 
medium (HIM), fir s t ident ified bv lField et all (j!969[ ) and 
iMcKee fe Ostrikerl (|1977|) . although in reality nm.ci is 
usually somewhat smaller than the chosen 1 atom cm~ 3 . 
Finally, in the fiducial model, all photons arc emitted in 
the line center, from the ICM in the center of the sphere. 

The covering factor f c is analogous to an optical depth 
of clouds intercepted by a sightline from the center of the 
galaxy and out. With a number density n c \ = N c \/V g& \ 
of clouds, the covering factor of the fiducial model is 
then f c = n cl r gal cr c i = |iV cl (r cl /r gal ) 2 ~ 2 (where 
a c \ is the cross section of a cloud). Since an aver- 
age sightline passing through a cloud traverses a dis- 
tance (d) = V c \/<j c \ = 4r c i/3 inside the cloud, the to- 
tal path traveled inside clouds for an average sightline is 
(din clouds) = fc(d) = iV cl r 3 /r 2 al . 

At the heart of the EW boosting mechanism lies the 
assumption of the Lya photons being shielded from the 
dust by neutral hydrogen while the non-scattered FUV 
continuum radiation is subject to the "full" dust ab- 
sorption. Hence, the relevant quantity for dust absorp- 
tion is the total average absorption optical depth (r a ) 
of the dust throughout the galaxy. The dust absorp- 
tion optical depth of a single cloud (ccnter-to-face) is 
T a c i = (1 — A)rd, c i) where rd, c i is the total (absorp- 
tion + scattering) dust optical depth, and A is the dust 
albedo. Hence, (r a ) ~ r a . cl (d in clouds) /r c i (this expression 
is only approximate, since photons scattering on dust 
grains may still be subject to absorption by another dust 
grain before it escapes the galaxy) . 

Dust grains are built from metals, and the density of 
dust is thus assumed to scale with metallicity and hydro- 
gen density. Since observationally the optical depth is 
not readily measured, we will instead refer to the metal- 
licity Z, which is usually more easily probed. Depend- 
ing on the actual extinction curve used, the constant 
of proportionality will differ. Furthermore, assuming 
that the density of dust scales linearly with the metallic- 
ity, the dust optical depth can be related to Z through 
t c i = A^HO'd(^) Z/Zq, where Cd is the dust cross section 
per hydrogen nucleus — specific for a given extinction 
curve — and Zg is the reference metallicity of that ex- 
tinction curve. For example, for SMC dust, at the Lya 
wavelength the cross section is <7d = 4 x 10~ 22 cm 2 with 
only a small wavelength depence across the line, and the 
SMC metallicity is Z ~ 0.25Z Q . 

Having explicated the basics of the model, we now pro- 
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Fig. 2. — EW boost 6 as a function of covering factor / c , for 
various values of cloud Hi densities njj : c j and metallicites Z. Lines 
with the same color share the same value of dust optical depth 
per cloud: blue, green, and red denote r a c i ~ 0.02, 0.2. and 2, 
respectively. For these simulations, the number of clouds and the 
galaxy radius was JV c i = 100-50 000 and r ga i = 2-10 kpc. Densities 
are given in cm -3 , and metallicites in terms of Zq. 

ceed to investigate the impact on the boost of varying the 
implied parameters. 

3.1. Cloud covering factor 

For the Neufeld scenario to be efficient, the sky, as seen 
from the center of th e galaxy, must be sufficiently cov- 
ered by clouds. IN91I asserts that the covering factor f c 
must be larger than unity. In fact, as soon as a single 
cloud is present, the 4tt average will have b > 1. Figure 
[2] shows, for various values of ram.ci and Z how the EW 
boost changes with f c . The boosts are seen to follow an 
approximate log-normal relation, where different combi- 
nations of riHi,ci and Z with equal dust optical depths lie 
roughly on the same line. At fixed Z, larger Hi densities 
result in larger boosts. The Lya photons are exposed to 
the same amount of extinction, but since (r a ) is larger, 
the continuum is reduced more. At fixed riHi.ci, a larger 
Z also results in a larger boost due to the higher (r a ), 
but since for each cloud interaction a Lya experiences a 
larger probability of being absorbed, combinations with 
high Z have smaller boosts than those with high 7iHi,ci- 
This is especially true at large values of f c , where the 
Lya photons interact with many clouds. 

3.2. Cloud hydrogen density 

To shield the Lya photon from the dust, the clouds 
must be highly optically thick in hydrogen. Since non- 
absorbed photons are usually reflected after only a hand- 
ful of scatterings, in principle the clouds could have an 
optical depth of only a few, although in that case, for 
realistic metallicites, the dust optical depth would be so 
low that the continuum passes unhindered through the 
clouds. 

Figure [3] shows how the magnitude of the boost de- 
pends on the Hi (column) density of the clouds. Note 
that for a given series of simulations, the metallicity has 
been held fixed, meaning that (r a ) increases with increas- 
ing iVHi.ci- Even for Solar metallicity, the EW boost is 
seen to set in only after Nui, c i surpasses approximately 
10 20 cm -2 , while, say, ten times lower metallicity re- 
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Fig. 3. — EW boost b as a function of neutral hydrogen column 
density Afjj : c i of the individual clouds, as measured from the center 
to the surface of a cloud, for different gas metallicitcs. The upper 
x axis gives the corresponding number densities in cm -3 in the 
case of r c i = 100 pc. Metallicities are given in terms of Zq . The 
simulations for which f c ~ 4 have been realized simply using JV c i = 
13 000. 
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Fig. 4. — EW boost 6 as a function of gas metallicity Z — which 
can be taken as a proxy for (r a ) — for different Hi densities njjj cl 
of the clouds. 

quires ten times higher JVhi, c i- The critical value for a 
noticeable boost is found where -/Vhi, c i and Z conspire to 
a dust optical depth of order unity. For very high values 
of -/Vhi,c1 the clouds become so optically thick in dust to 
continuum such that the calculated EWs become rather 
noisy. 
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Fig. 5. — Lower x axis and red lines: EW boost b as a function 
of cloud radius r c i, maintaining constant covering factor f c (solid) 
and a constant total Hi mass (dashed). Upper x axis and green 
line: EW boost as a function of cloud size distribution power-law 
index /3, for minimum and maximum cloud radii r m i n = 20 pc and 
r maJ = 200 pc, at constant f c . 

the adopted cloud size of 100 pc, this corresponds to the 
product nui Z/ Zq ~ 1, which is indeed seen in Fig. 2] 

3.4. Cloud size distribution 

In the previous simulations the c loud radii wer e held 
constant at r c \ = 100 pc. As IN91I proposed and IHO06I 
confirmed through their numerical model, the covering 
factor rather than the actual shape and size of the clouds 
is the important parameter controlling the RT of the 
Lya photons. For the continuum, however, doubling 
the cloud radius at constant / c doubles the dust optical 
depth, thus increasing the boost. In Fig. [5] we investigate 
the dependency of the boost on cloud sizes, as well as on 
cloud size distributions P{r c \) oc r~f . The dependency 
on slope (3 is seen to be small. In the simulations with 
varying cloud radii, the slope wa s held constant at j3 = 2. 
This value was chosen to follow IHO 061 . 

Increasing r c \ or decreasing /? at constant / c both corre- 
sponds to in increasing the (tj) and hence b, although it 
is seen that changes in (3 affects b only minutely. On the 
other hand, changing the cloud sizes but maintaining a 
constant total Hi mass implies a constant (Td), and hence 
a constant continuum escape fraction. Consequently, the 
boost does not change for r c \ > 100 pc, but at lower r c i 
the number of times that the Lya photon interact with 
clouds (-/Vo) becomes so large that the Lya escape frac- 
tion decreases considerably, reducing the boost. 



3.3. Dust contents in the clouds 

The density of dust scales both with metallicity and Hi 
density. Increasing the gas density in the clouds will in- 
crease the dust density accordingly. This does not change 
the surface of a cloud from the point of view of a Lya 
photon, but the continuum photons experiences a larger 
total optical depth of dust, and thus the boost increases. 
For a given Hi density, increasing the metallicity results 
in a higher boost. These relations are shown in Fig. @] As 
expected, larger Hi densities result in larger boosts, but 
for sufficiently low metallicitcs the boost will disappear. 
Again, the critical value of Z is the metallicity required 
for the optical depth of dust to be of order unity. For 



3.5. Cloud velocity dispersion 

The very concept of the boost hinges on the fact that 
neutral gas shields photons close to the line center from 
the dust. But if a given cloud as a whole has a non- 
vanishing bulk velocity, the entire spectrum is Dopplcr 
shifted in the reference frame of the cloud, and for suffi- 
ciently large velocities, all line photons suddenly are no 
longer in resonance. In general, Hi regions may be ex- 
pected to exhibit such macroscopic motions, described 
by a cloud velocity dispersion o~v,ci- 

Figure [5] demonstrates that this is indeed the case; for 
°V.ci ^ 100 kms -1 , the fiducial model is seen to be in- 
capable of boosting the EW. Moreover, whereas in the 
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Fig. 6. — EW boost as a function of cloud velocity dispersion 
<j Vc i, for zero-width (solid) and broadened (dashed) emission, in 
galaxies of f c ~ 2 (green) and 4 (red). 
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Fig. 7. — EW boost as a function of cloud gas temperature T c i 
(including also small scale turbulence) for a covering factor f c ~ 2 
(green) and / c ~ 4 (red). 

case of ovci = a higher covering factor results in a 
larger boost, introducing random cloud motion reduces 
the boost even faster; for f c twice that of the fiducial 
model, the boost vanishes already for o~y c \ > 30 kms . 



3.6. Cloud gas temperature 

Increasing the gas temperatur^l T c \ in the clouds 
means less atoms with the right velocity for scattering 
photons exactly at the line center, in return for more 
atoms available for scattering off center photons. Pho- 
tons in the wings of the line, however, do not care about 
the gas temperature, since the profile here is given by 
natural broadening. Even so, the difference in the effec- 
tive cross section of the atoms is not very large. The 
consequence is that, for increasing T c \, the photons pen- 
etrate slightly deeper into the cloud, but not enough to 
affect the boost, as seen in Fig. [Jj The only difference is 
that the line profile of the escaping radation is somewhat 
broadened. 

3 The parameter "temperature" covers also sub-grid turbulent 
motion, since this effectively broadens the line profile by a Gaus- 



FlG. 8. — EW boost as a function of terminal outflow expansion 
velocity Voutfor covering factors / c ~ 2 (green) and ~ 4 red. 

3.7. Galactic outflows and infall 

Gas elements also move in large-s cale, collective mo- 
tions , se en e.g. in ga lac tic ou tflo ws dJohnson fc Axfordl 
19711: IRiroke et al l 1200% IShaolev et all Effll 
Veilleux et al.l I2005D . For the same reason as in 
the previous section, such motions should act such as to 
diminish t he boost. This stands in contrast to what is 
expected (jKunth et al.l 119981: IDiikstra fc Wvithd 12010ft 
for a more homogeneous shell of gas being expelled from 
the galaxy, surrounding a central source and covering 
the full sky. In that case, the Lya photons that would 
otherwise have to scatter their way through the shell, 
being very vulnerable to dust absorption, may be shifted 
away from the line center and escape through the shell 
with minimal absorption. 

Figure [8] shows how the EW boost depends on the 
expansion velocity V ou t- The wind speed is a func- 
tion of distance r from the center and is modelled as- 
sumed th at the gas elements receive an acceleratio n 
a cx r~ a (jSteidel et alJl20l5t IDiikstra fc Kramerll2012l ). 
The speed thus increases from at r — to the terminal 
velocity V^ u tat r = r ga i. Note, however, that the exact 
wind profile is not a major determinant of the resulting 
boost. 

The curves in Fig. [8] resemble those in Fig.[6l although 
they are somewhat more shallow. The reason is that a 
photon bouncing off of a cloud in a galactic wind will, 
in general, upon its next encounter with a cloud, meet a 
cloud with roughly the same velocity. In the case of ran- 
dom motions, the relative velocity of the next cloud may 
be very high, increasing the probability of the photon 
being absorbed. 

For reasons of symmetry, in terms of b the same results 
are obtained if the velocities are in yerted, as would be th e 
case for galactic gas accretion (e.g. lDiikstra et al"1l2006 bl ) . 
The spectra, however, would be reflected about the line 
center. 

3.8. ICM hydrogen density 

Thus far, we have assumed an ICM completely void of 
both dust and neutral gas. Due to the high temperature 
of the ICM, the hydrogen is expected to be quite highly 
ionized. Similarly, due to the high temperature, as well 
as the stronger ionizing UV radiation field, dust may be 
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Fig. 9. — EW boost b as a function of Hi density JVjj : icm m * ne 
ICM. Results are show in red, green, and blue for cloud number 
densities njjj cl = 0.1, 1, and 10 cm -3 , respectively, in combination 
with cloud metallicites Z = 0.1 ZQ(dashed) and 1 Zq (solid). 

expected to be partly depleted. However, even a single 
scattering in the ICM may be fatal to a Lya photon: Ex- 
actly because of the high temperature, a scattering event 
in the ICM is likely to occur on a high-velocity atom 
moving more or less perpendicular to the path of the 
photon, such that the velocity relative to the trajectory 
of the photon is small. Unless the photon is scattered 
close to (or opposite to) the same direction, it will be 
highly Dopplcr-shifted, such that next time it encounters 
a cloud, it will penetrate deeply into the dusty medium, 
with a higher probability of being absorbed. This effect 
is investigated further in Sec. 13.101 

For our fiducial model's radius of r ga i = 5 kpc, an 
optical depth of order unity in the ICM for a line center 
photon is reached for an Hi density Uhi ~ cm -3 . 
Due to the resonance nature of the scattering, somewhat 
higher densities are allowed as soon as the photon has 
diffused a few Dopplcr width from the line center. This 
effect is visible if Fig. [9j where the boost is shown as a 
function of ICM Hi density for different cloud densities 
and metallicites. 

3.9. ICM dust density 

The neglect of dust in the ICM may be justified by the 
fact that dust grains will tend to get destroyed in hot 
and ionized media. Nevertheless, dust does exist in Hn 
regions, albeit generally at lower densities. For an in- 
creasing covering factor the photons' paths through the 
ICM is increased as they walk randomly out through the 
galaxy. The mean free path between each cloud interac- 
tion is l// c , while the average number of cloud interac- 
tions is Nq — f c + |/ c (see App. IA.4]) . Thus, the av- 
erage path length for a photon out of our fiducial model 
01 r gai = 5 kpc is (dtot) — 14 kpc. For a density of 0.01 
cm~ 3 , Z = Zq, and no dust destruction, the total optical 
depth of dust actually reaches (ticm) ~ 0.2. For larger- 
densities r ga i and / c , (ticm) could exceed unity, while if 
part of the dust is destroyed, it may be neglected alto- 
gether. 

Figure [TU] shows the effect of dust in the ICM for var- 
ious densities, galaxy radii, and dust depiction factors. 
The important quantity is the number density of dust, 



Fig. 10. — EW boost b as a function of dust density in the inter- 
cloud medium, parametrized through the amount of metals in the 
form of dust, (Z/ Zq) x (^ icm/€c1 o)i where the last term is the 
dust-to-metal ratio in terms of the values typically found in the 
predominantly neutral environments from where the dust extinc- 
tion curve is taken (here the SMC). Results are shown for galaxies 
°f T 'gal = 5 (green) and 10 kpc (red), in combination with total 
(neutral + ionized) hydrogen density «h,icm = 10~ 3 (dashed) 
and tihjcm = 10 -2 cm -3 (solid). 

which is proportional to the product of total density, 
metallicity, and dust-to-metal ratio £d! that is, a given 
dust density can be realized in several ways. While in 
predominantly neu tral regions, £ d is rather constant (see 
discussion in Sec. 14.41 and Sec. 14. 8[) . Accordingly, the 
boost is shown as a function of the product of Z and £d , 
for different total hydrogen densities. Indeed, ICM dust 
is seen to be important only for rather high ICM densi- 
ties, large galaxies, and low dust depletion factors. Thus 
we see that the most crucial factor of the Neufeld mecha- 
nism is not, as it is often laid out, that the Lya photons 
are confined to a dustless medium but rather that are 
confined to a medium of low density but possibly with 
the same dust-to-gas ratio, and high ionization fraction. 

3.10. ICM gas temperature 

As mentioned in Sec. 13.81 although scatterings in the 
ICM evidently are rare a single such event may easily 
result in the subsequent absorption of the photon. Fig- 
ure QT] explores this effect, but it is seen that only for 
rather high values of tihiJCM, the boost decreases sig- 
nificantly. For low densities, scatterings are so rare that 
even though some photons are absorbed, it does not af- 
fect the net result much. For very high temperatures, the 
boost increases again, since the number of atoms avail- 
able with the right velocity to scatter a photon becomes 
too small. 

3.11. Location of emission 

Al though star formatio n tends to be centrally located 
(e.g. lFruchter et aT1l2006|) . obviously the photons are not 
emitted from a point source in reality. However, as 
long as the initial emission direction is isotropically dis- 
tributed, photons emitted in the outskirts of the galaxy 
may be emitted both in an outward direction, being sub- 
ject to a lower covering factor and thus a smaller boost, 
or toward the center, resulting in a higher covering factor 
and higher boost. As the number of cloud interactions 
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Fig. 11. — EW boost as a function of temperature Tjcm in the 
ICM, for neutral hydrogen densities njj : ICM = 10 (dotted), 

n Hi ICM = 10 -6 (dashed), and njjj ICM = 10~ 4 (solid) cm -3 . 
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Fig. 12. — The effect of an extended emission, as given by an 
exponentially decreasing emission site PDF of the form P(r) oc 
e~ r / H , where r is the distance from the center of the galaxy, and 
H is the scale length. The rightmost point shows the boost for 
H — > oo, i.e. for a completely homogeneous distribution of sources. 

scale non-linearly with / c , however, an extended emission 
profile does result in a slightly lower boost, although the 
effect is miniscule, as seen in Fig. [T2] 

While the extended emission relaxes the notion of a 
central point source, all photons were still emitted from 
the ICM, as in the case of the central source in all pre- 
vious sections. However, as the Lya radiation is as- 
sumed to originate from gas surrounding young stars, 
which form from gas clouds that have recently cooled 
sufficiently to initiate star formation, it is probably more 
realistic to have at least a fraction of the photons origi- 
nating from within the clouds. In Fig. [13] we investigate 
how the EW boost changes for an increasing correlation 
of the photon sources with the clouds. As expected, tak- 
ing this effect into account diminishes the boost, as the 
path out of a cloud is longer for a Lya photon than for 
a continuum photon. 

3.12. Intrinsic line profile 

In the fiducial model, all Lya photons are emitted ex- 
actly in the line center, i.e. as a delta function. In reality, 



Fig. 13. — EW boost b as a function of probability P r \ for a 
photon to be emitted from inside a cloud. An emission PDF with a 
scale length of H = 1 kpc has been adopted. For P c \ = 0, the model 
reduces to the fiducial model, except for the somewhat exte nded 
emission which has only a modest influence on b (cf. Fig, I12[l . For 
P c l = 1. all photons are emitted from the interior of clouds. Results 
are show for all combinations of / c ~ 2, 4, njj : cl /cm~ 3 = 1, 10, and 
Z/Zq = 0.1, 1. 
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Fig. 14. — EW boost b as a function of line width cr em for covering 
factors f c ~ 2 (green) and 4 (red). 

the intrinsic line will have a finite width, which is con- 
volution between the natural, Lorentzian line profile and 
the thermal, Gaussian distribution of the atom velocities. 
Moreover, macroscopic motion of the gas parcels emitting 
the photon will cause an effective broadening of the line. 
As the temperature of the emitting gas is of the order 
10 4 K, corresponding to a line width of ~ 10 kms -1 , 
the macroscopic motion may dominate over the thermal 
(see Sec. 14. ip . From non-resonant lines such as Ha, we 
know that the intrinsic line width may easily reach tens 
of kms -1 , i.e. several Doppler widths. As photons born 
this far from the line center are not efficiently shielded 
from the dust, they are more prone to dust absorption. 

The effect of a broadened emission line profile is shown 
in Fig. [H 

4. TOWARD A REALISTIC MODEL 

Whereas the fiducial model in the previous section 
was not chosen to be particular realistic, but rather to 
roughly resemble the original Neufeld scenario, we now 
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proceed to investigate what values of the physical param- 
eters we in fact expect to meet in real LAEs. In building 
a "realistic" model, it should be kept in mind that sev- 
eral of the input parameters are mutually dependent. For 
instance, under the assumption of approximate pressure 
equilibrium, the temperatures and densities are related 
as rctot.cfTci ~ «tot,iCM7icM, where subscripts "tot" refer 
to the total number density of all particles. The Hi den- 
sities are then given by the ionization fractions which, 
in turn, depend on not only temperature, but also the 
chemical composition, i.e. on metallicity. Note, however, 
that due to very different cooling time scales of the WNM 
and HIM, the different phases of the ISM need not be in 
thermal equilibrium. Furthermore, the dependencies of 
some of the parameters may actually serve such as to 
counteract a boost. For instance, a larger galaxy in gen- 
eral entail a larger / c , which may increase b, but since 
the larger mass also typically implies a larger oy, c i, b is 
quickly reduced again. Similarly, although a increased 
metallicity tends to increase b, metals also donate free 
electrons, which aids to increase the neutral fraction of 
hydrogen in the ICM, thus decreasing b. 

The approach taken here is to identify for each param- 
eter the minimum and the maximum values conceived to 
be encountered in LAEs, as well as more typical ranges. 
Subsequently a series of RT simulations is run with the 
physical parameters taking random values in the chosen 
ranges, and finally regions of high and low likelihood are 
identified in the parameter space, exploring whether any 
of these models lead to a significant boost. 

LAEs span a wide range of physical characteristics, and 
in many aspects it does not make sense to speak about 
a "typical LAE" . It is often stated that LAEs are quite 
young, fairly small, relatively highly star-forming galax- 
ies of rather low metallicity and dust contents. However, 
one must bear in mind that many of these aspects are a 
consequence of the selection criterion picking out objects 
at high redshifts, i.e. at early epochs. In the following 
sections, values and ranges of the relevant parameters 
argued to be more or less representative of LAEs are 
summarized from the literature. An overview of the val- 
ues is found in Tab. [1] on pagefT^l 

4.1. Cloud hydrogen density and temperature 

Gas does not cool equally well at all temperatures. 
Rather the cooling function exhibit several plateaus at 
which the gas temperature tends to settle. This is part 
of the reason for the different phases of the ISM; for 
gas which is in rough pressure equilibrium, the densities 
of the phases are set by their temperature, and thus, 
although variations do exist, they are not expected to 
span many orders of magnitude within a given phase. 
In the WNM, to which the clouds in the model corre- 
spond, characteristic hydrogen densities are 0.2-0.5 (e.g . 
ICarilli et alll^lFerrierel2001tlGloeckler fc GeissJ200l 
— typically probed by Hi 21 cm emission — with a neu- 
tral fraction close to unity. At low (high) pressure, njji 
may reach values of ~ 0.02 (3) cm- 3 (|Wolfire et al.|[2003h . 
If densities rise much above ~ 1 particle cm~ 3 , cooling 
becomes so efficient that the gas will tend to contract 
further to a cold neutral medium (CNM) and molecular 
clouds, thus quickly decreasing f c . The temperature of 
the WNM is of the order 10 4 K, ranging from 8 000 K to 



12 000. Since small scale turbulence effectively broadens 
the line in the same way as thermal motion, this pro- 
cess can be taken into account using a higher tempera- 
ture. Turbulence is usua lly of the order of (a few times) 
the speed of sound (e.g. IGaensler et al.|[20lTl ) which, in 
turn, is the same order as the thermal motion. The ef- 
fect can thus be simulated by letting T = 2-3 x 10 4 K. 
Note also that much of the WNM may b e in a thermally 
unstable phase (|Heiles fe Trolandl |2003|) . reaching tem- 
peratures below 5 000 K. 

4.2. ICM hydrogen density and temperature 

Surrounding the neutral gas is the HIM, heated and 
ionized to a large degree by supernova shock waves 
sweeping through the ISM. Expanding due to the larger 
pressure, and thus diluting, the cooling timescales are so 
long that temperatu res of 5 x 10 5 -5 x 10 6 K and higher 
is easily maintained (jBrinks et al.ll20dot iTiillmann et al.l 
120061 120081 ). Densities in the range 10~ 3 to 10~ 2 ar e 
typical here (jDopita fc Sutherland! [2001 lFerrierel[200l . 
but since the fraction of neutral hydrogen at thes e tem- 
perat ures is of the order x ht.tom ~ 10~ 8 -10~ 5 (|Housd 
11964 ISutherland fc Dopitolll993l ). scattering in the ICM 
is very rare. These pr operties are usually pro bed through 
X-ray emission (e.g. ITiillmann et al.l 120061). or absorp - 
tion lines of highly ionized metals (e.g. iPark eFaDEOOg). 
Note that the shock waves push the phases out of equi- 
librium and may in fact dominat e the state of the ISM 
(|de Avillez fc BreitschweTdt|[2005l ). 

4.3. Covering factor 

A typical covering factor is difficult to establish, but a 
hint is provided by spectra of s ources lying behind inter - 
vening galaxies. For instance, iNoterdaeme et al] (|2012D 
fit low-ionization absorption lines such as Sin — orig- 
inating from cool clouds — in a quasar spectrum, the 
line of sight toward which passes through a DLA lying 
at z = 2.2, with 5-8 components. The impact parameter 
of the sightline with respect to the central emission of the 
DLA galaxy counterpart is pretty small (< 1 kpc), and 
thus in this case the covering factor is roughly f c = 2.5-4. 

In the Milky Way (MW), the numb er of neutral clouds 
per kpc in t he plane is roug hly 4-8 dKnudd 119791 li~98l 
lFrancoll2012D . Although this implies a high covering fac- 
tor along a sightline in the plane, perpendicular hereto 
the covering factor is only of order unity. Since the disk 
has formed from a collapsed sphere, one may expect a 
smaller number of clouds per kpc for a spherical system. 
For instance, the scale height of the MW's WNM is ap- 
proximately h = 300 pc an d its radius is i?MW ~ 15 kpc 
(jBinnev fc Merrifieldlll998l ); distributing the clouds over 
a sphere of the same radius reduces the crowdedness of 
clouds by a factor of ~ Ryiw/h = 50. 

4.4. Cloud metallicity and dust density 

Generally lying at high redshifts where metals have 
had less time to build up, LAEs tend to have lower 
metallicities, and hence dust densities, than local galax- 
ies. Measured values of Z range from metallicities as low 
as 1§~ 2 Zq to values similar to local values, and in 
some cases even super solar metallicities have been found 
at high redshifts (e.g. lAcauaviva et a l. 2012). In a sur- 
vey of ~ 10 3 LAEs at z = 2.2, iNakaiima et al.l (|2012[ ) 
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found a lower limit for the average mctallicity of a LAE 
of Z = 0.09Z Q . Even low-redsh ift LAEs tend to have 
modest metallicities (Z ~ 0.1-1 iGiavalisco et al.l IT9961 : 
IQstlin et al.ll2009h . 

Correspondingly, measured col or excesses tend to be 
mo dest. For example, at z ~ 3 iGronwall et al.l (|2007l ) 
and lPirzkal et al.l (|2007T > found E(B - VYs lying approx - 
imately in the range 0.01 to 0.1. iVerhamme eTaLl ((2008) 
found similar values by fit ting Lyct line profiles using 
the Lya RT code MCLya (jVerhamme et al.l 120061 ) . At 
z = 2.2. lHaves et al.l (|2010T l found values spanning all the 
way from to ~ 0.3. For comparison, our fiducial model 
has E{B -V) = 0.05. 

In our model we assume that the dust-to-metal ratio 
£d,ci in the clouds is equal to that of our reference extinc- 
tion curve, such that Td scales linearly with Z . Observa- 
tions of £d,ci at high redshifts a r e sparse, but tend to be 
similar (e.g. IPettini et al.l 119971: iSavaglio et al.ll2003D or 
slightly lower. Various analytical and numerical calcu- 
lations predict a metallicity-dependent evolution of the 
dust-to-metal ratio, such that £d,ci reaches present-day 
values only after timescal es of 10-10 Myr dGall et al. 

mor even several Gyr (|Inouell20'03l see also IMattsson 
). Moreover, measurements in the local but very 
low-metallicity galaxy I Zw 18 indicate that £d,ci is lower 
in such galaxies. Since we have used the extinction law 
of the SMC, which is also a low-metallicity galaxy with 
a rather young stellar population, we expect that its ex- 
tinction will not be substantially different from that of 
high-rcdshift galaxies. 

4.5. Cloud velocity dispersion 

From the virial theorem, the components of a galaxy 
of mass M and radius r wil l have ch aracteristic veloc- 
ities of the order ov, c i = \J GM/Cr, where the factor 
C depends on the geometry and the actual mass dis- 
tribution, as well as on whether the system is rota tion- 
or dispersion-dominated (jBinnev &: Trema inc 20081) . For 
dispersion-dominated galaxies, M refers to the total, 
dynamical mass of the system, and C ~ 6.7 for vari- 
ous g alactic mass distributions (jForster Schreiber et al.l 
I2009D . For rotation-dominated galaxies, the appropriate 
mass is the mass enclosed within r, and C ~ 2.25, again 
averaged over vario us galactic mass distribution models 
(jEpinat et al.ll2009f ). 

Measured values of the Hi velocity dispersion, or of the 
stellar velocity dispersion which arguably reflects the ve- 
locity field of the gaseous compone nts, range from ~ 10- 
30 km s" 1 for dwarf gala xies (e.g. iPeterson k, Caldwelll 
119931 : Ivan Zee et al . 1998) to severa l hundred kms" 1 for 
large ellipticals (e.g. lMcElroy||1995|) . Note that observa- 
tionally what is measured is the dispersion along a line of 
sight, or an average of many lines of sight, and hence cor- 
respond to the dispersion in one dimension. If there is no 
preferred direction of motion, the three-dimensional ve- 
locity dispersion, which is what is referred to in Sec. 13.51 
is a factor s/3 times higher. 

Our fiducial model has a total Hi mass of 6.7 x 10 8 Mq. 
For small galaxies, t he ratio of Hi to the total, dynamical 
mass is roug hly 0.1 (ISkillman et al.lfl987t Ivan Zee et al.l 
119971 fT998h . Thus, the velocity dispersion should be ap- 
proximately 50 kms" 1 . 

At high redshift, most galaxies may not yet have had 



time to form a disk, and may thus be expected to be 
dispersion-dominated. For galaxies which eventually set- 
tle i nto a disk, the veloc ity dispersion tends to decrease 
(e.g. lThomas et all 20121). For instance, fo r intermediate- 
to late- type spirals. iBershadv et al.1 (|2011l ) found that the 
central vertical velocity dispersion is ~ 1/4 of the maxi- 
mum rotation speed. 

The preceding reasoning dealt with the overall veloc- 
ity dispersion of the galaxy's components. For relaxed 
disks, these could be correlated in phase space, such 
that locally to a given photon's emission site the first 
few clouds encountered will have a lower relative veloc- 
ity. With a lower limit of 5 kms , however, we believe 
that a realistic threshold has been met; for instance, in 
the MW the dispersions in peculiar motion of local stars 
with respect to the loc al standard of rest is ~ 20 kms -1 
(jSchonrich et al.ll2010h . 

The velocity dispersion of the gaseous component 
might be expected to be smaller than the stellar veloc- 
ity dispersion, however, but in The Hi Nearby Galaxy 
Survey (THINGS), the galaxies, which are almost exclu- 
sively disk galaxies, exhibit velocity dispersions between 
10.1 and 24.3 kms -1 , with a mean of 16.8 ± 4.3 kms -1 
(jlaniamasimanana et al.ll2012f ) . 

4.6. Galactic outflow velocity 

Although not ubiquitous, galactic winds seem to be 
rather common in high-rcdshift galaxies. The outflows 
are produced by the kinetic and/or thermal feedback 
of massive stars and supernovae on the ISM. Radia- 
tion pressure and heating engender expanding bubbles 
of primarily ionized gas that sweep up interstellar ma- 
terial, eventually escaping the galaxy (or possibly re- 
enterin g the system a fter having reached several virial 
radii) ( Heckmanl 120021 ) . Outflows are usually detected 
through low-ionization absorption lines such as Mgn and 
Fell which appear blueshifted with respect to the sys- 
temic |_v^kodty^_Tjykalvaluesa^ 100 kms -1 
(e.g. IPettini et alt 120011 : iRubin et all 120111 ). but range 
all the way from and up to ~ 1-2000 kms -1 (e.g. 
IShaplev et all[2003h . 

Various models fo r gala ctic outflows have been put 
forward: lOrsi et all (|2012f ) equate Voutwith the circu- 
lar velocity, which is of the same order as oy. r i, dis- 
cussed in th e previ ous section. iBertone et al.l (|2005l ) and 
iGarel etldl (|2012f ) apply a weak dependency on SFR, 
with Kut cx SFR 0145 , where the SFR is given is M Q 
yr _1 , and the constant of proportionality is ~ 300-1000 
kms -1 . 

In the context of LAEs, the effect of an outflow is 
to diminish the blue peak of the otherwise symmetric, 
double-peaked line profile. This is noticable already ~ 10 
kms" 1 , and by ~ 10 kms -1 , the blue p eak may be 
erased altogether fe.g. lDiikstra et ai1l2006al ). Most high- 
rcdshift LAE profile observations have resulted in the 
well-known, asymmetric red peak only. However, using 
sufficiently high resolution, R > 1-2000, it seems that a 
significant fraction (20-50%) of LAEs actually show at 
least signatures of a blue peak, indicating only modest 
outflow velocities dVenemans et al.l 120051: iTapken et al.l 
l2007llKulas et al.l 120 12HYamada et al. f l2012D . 

This hypothes is is in qualitativ e conc ord with the re- 
cent findings of lHashimoto et ail (|2012l ). who interpret 
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the smaller offset of LAE Lya lines with respect to the 
systemic velocity, compared with those of LBGs, as be- 
ing due to smaller outflow velocities. Note, however, that 
the line offset to a large degree depends on the column 
density, which due to the generally much higher gas mass 
of LBGs naturally will be larger. 

For the reasons given above, and since a large outflow 
velocity was found to quicky destroy the boost anyway, 
we will restrain ourselves to investigating rather small 
values of V ou t, from to 100 kms" 1 . 

4.7. Emission sites 

Since evidently the extend of the photon-emitting re- 
gion is not of major importance (see Sec. 13. lip . for sim- 
plicity we will confine our grid of models to one value 
only, viz. H cm = 1 kpc. This value is roughly equal to 
the UV half-l ight radius of LAEs a t redshifts 2 < z < 6 
(see Fig. 2 of lMalhotra et al.ll2012t) . 

The typical environment from which photons are emit- 
ted is more crucial. Originating in the Hi/Hii-boundarics 
surrounding massive stars, the concept of emitting the 
Lya photons from the ICM is rather dubious. Mas- 
sive stars are predo minantly born in gian t molecular 
clouds (GMCs; e.g. IGarav fc Lizanol |l999) Since such 
stars are short-lived, with lifetimes of but a few Myr, 
they are not expected to travel very far from the dense 
clouds from which they were born. This favors a high 
value of P c i, the probability of being emitted from a 
cloud rather than from the ICM. On the other hand 
the intense UV radiation carves out a Stromgren sphere, 
whic h reaches a size of the order of 10(s) pc in a few 
Mvr dHosokawa fc Inutsukall2006l:lGendelev fe Krumholzl 
I2012T ) . If they are matter-bounded rather than radiation- 
bounded, the Lya photons start to stream freely out into 
the surroun ding, l o wer-d ensity medium. 

Although llsraell (119781 ) found that most massive stars 
are located i n the outskirts of GMCs (in the MW), 
IWaller et al.l (fl987T) argue that they tend to be more cen- 
trally located, but that ~ 30% of the Hii regions used to 
probe O and B stars in reality are just parcels of ionized 
gas that was once in molecular for m near the surface of its 
host c loud. On the other hand, iGendelev fc Krumholzl 
reason that internal turbulence in the GMCs cre- 
ates filamentary structure so that any star has a high 
probability of being born near the surface. 

Regardless of the massive stars being born deep in the 
GMCs, or being born closer to the edge so that blisters 
allow a high escape fraction of ionizing and hence Lya 
photons, the GMCs themselves are usually embedded in 
the WNM. The density of the WNM being several orders 
of magnitude lower than the GMCs, a Stromgren sphere 
that breaks out of a GMC may be able to grow faster, 
faciliting escape into the ICM. 

Massive star do also exist in exposed cluster, such 
as in the Pleiades. Although the stars necessarily 
are born in dense clouds, some stellar groups man- 
age to ionize and blow away the neutral gas. Obser- 
vations of this at high redsh ifts are difficu l t, bu t in 
the nearest 2 kpc of the Sun, lLada fc Ladal (|2003l ) es- 
timate the (lower limit) birthrate of embedded clus- 
ters to be 2-4 Myr -1 kpc -2 , which is more than an 
order of magnitude highe r than that of open c l usters 
(0.25-0.45 Myr" 1 kpc 2 lElmegreen fc Clemens! U985I: 
iBattinelli fc Capuzzo-Dolcettalll991l) . 



In a study of 45 GMCs in M33, llmara et al.l (|20T1 
found 29 to be spatially and kinematically coincident 
with a local peak in atomic gas, 13 to be kinematically 
coincident, but located near the edge or on a filament 
between two peaks, and the remaining three not to be 
associated with any high-densit y atomic gas. In t he re - 
gion NGC 602/N90 in the SMC. lGouhermis et al.l ((20121) 
found an "unusually large fraction" of 60% of the pre- 
main sequence stars to be clustered, while the rest are 
diffusely distributed in the intercluster Rl OR. 

Additionally, in the MW 10-20% of all O stars are 
found in ultracompact Hii regions, still em bedded in their 
natal molecular cloud (|Chur chwcll 199(3). 

From these considerations, we consider P c i = 0.1 a 
lower value, 0.2-0.5 to be a more or less realistic value, 
and 0.9 to be a highest value. 

4.8. ICM dust density 

Since the various ways of destroying dust tend to corre- 
late with processes that also ionize gas (collision, sputter- 
ing, sublimation, and evaporation), the dust-to-metal ra- 
tio £dj an d hence the dust-to-gas ratio, in ionized gas may 
be expected to be lower than in neutral regions. How- 
ever, most extinction curves are obtained from sightlincs 
crossing several phases of the ISM, although typically 
only the correlation of the extinction with the neutral 
hydrogen is probed. Hence, the variation of £d with xhi 
is not well-constrained. Whereas £d hi neutral gas seems 
to be more or less universal over both different phases, 
metallicities, and redshifts, the observations that do ex- 
ist generally indicate a lower £d in regions of ionized gas, 
albeit with large variation s. On the basis of a discus- 
sion of various such regions. ILaursen et al.l (|2009bl ) argue 
(Sec. 2.1.1) that £d in these regions ranges between a fac- 
tor of 10 -4 lower than in the neutral phases, to values 
similar to these, but with typical values of 10~ 2 times 
lower. This is an average over several different types 
of Hii regions, in particular also the very dust-depleted 
IGM; in the HIM £d = 10 -1 may be more representative. 
A sim ilar discussion was given recently by iPaladini et al.1 
(|2012t Sec. 5), where also it is noted that the mere ra- 
diation pressure may be able to clear large volumes of 
dust. 

4.9. Emission line width 

The intrinsic Lya line width is given by a convolu- 
tion of the natural, thermal, and turbulent broadening. 
As discussed in Sec. 14.11 turbulence may dominate over 
thermal motions. Since scattering changes the line shape 
in a highly non-trivial way, Lya line shape observations 
do not reveal the intrinsic shape. Instead the width of 
non-resonant lines such as Ha, which are produced in the 
same locations as Lya, arguably can be used as a probe. 
Since the emitting regions also exhibit larger-scale mo- 
tions, however, when observing the line integrated over 
the full galaxy it will usually be much broader. 

When galaxies are well resolved, the internal motion of 
cloud are more e asily disentangled. In nearby galaxies, 
(jYang et al.lll994l ) found Ha line widths (in terms of stan- 
dard deviations) of 20-30 kms -1 . At higher redshifts, 
using CO lines where large-sc ale motions have been re- 
moved, [SwinbanF^F^L| ([201 ID determine the internal ve- 
locity dispersion of the components of a starburst galaxy 
at z = 2.3 to be 45-85 kms -1 . 



12 



Laursen et al. 



TABLE 1 

Summary of parameter values 



Parameter Full range 



"Typical" 



"Unusual" Sec. 



™Hi,d 

T A 

Z 

°V,cl 

fc 

'cl 

p 

"■Hl.ICM 

Vout 
Hem 
Pel 
C cm 

r ga i 



0.03-3 

5xl0 3 -3xl0 4 

0.03-2 

5-100 

0.8-8 

0.03-0.2 

1-2.5 

10 _12 -10 -6 

3xl0 5 -5xl0 7 

0-100 

1 

0-1 
5-100 
5, 10 



0.2-0.5 
8 000-20 000 
0.05-0.3 
30-50 
Full range 
Full range 
Full range 
10-iO-io-T 

1 B.5_ 10 6.S 

10-50 
1 

0.2-0.5 
10-85 
5, 10 



0.1-1 
Full range 
0.03-1 
10-80 
Full range 
Full range 
Full range 
Full range 
Full range 
Full range 
1 

0.1-0.9 
Full range 
5, 10 



1.1 



no 



go 



1.2 



■1.2 

TU 

4 
1.9 



Note. — Densities are given in cm , metallicity in terms of 
the Solar value, velocities and line widths in krns" 1 , tempera- 
tures in Kelvin, distances in kpc. 

4.10. Cloud size distribution 

As mentioned previously, the covering factor rather 
than the actual shapes and sizes of the clouds are the 
crucial factor in determining the boost. However, for 
P c i > 0, the size of the clouds become important, as 
smaller clouds implies an easier escape from the first 
cloud. On the other hand, at a given covering fac- 
tor a larger cloud size implies more attenuation of the 
continuum. In Sec. 13.41 we us ed = 2 when keeping 
r c i constant in order to follow IHO06I. Measured slopes 
tend to be more shallow ~ 1.6 (|Dickev fc Garwood! 
ll989HWilliams fc McKell997l ). Obviously, in reality the 
WNM does not consist of spherical clouds, so in order for 
the conclusions not to depend too much on the chosen 
cloud sizes, we will study a broad range of cloud sizes 
and distributions. Approximately half of the simulations 
will be run with fixed cloud size, both r c \ = 30 pc and 
r c i = 100 pc, while the other half is run with r c \ = 20-200 
pc and /? = 1-2.5. This cloud size range is rou ghly con- 
sistent with that of the LMC (|Kim et al.ll2003l) . Smaller 
cloud sizes probably requires further cooling, while larger 
clouds will get torn apart by la rge scale moti on, at least 
in rotation-dominated objects (|Newtonlll980f h 

5. RESULTS 

Having settled upon more or less realistic ranges of the 
parameters important for the RT, summarized in Tab.[TJ 
a large number (4 x 10 4 ) of simulations is carried out in 
the many-dimensional parameter space spanned by 
the parameters. Whereas in all previous simu- 

lations 10 photons were used to ensure convergence, in 
order to be able to adequately sample & , we use only 
10 3 photons per model. In most simulations, the result- 
ing boost is accurate to within roughly 10%; in simula- 
tions where the continuum escape fraction is very low, 
usually implying a large boost, this tends to overesti- 
mate the boost. One hundred of such simulations were 
resimulated to convergence; of these, five were found to 
have underestimated b by 20-30%, while the rest had 
overestimated b by ~ 0-300%. However, as will be evi- 
dent below, all of these models correspond to physically 
extremely unlikely models. 



The results of the simulations are displayed in Fig. [T31 
where the boosts of all models are shown as a function of 
total, average column density (Ahi) of neutral hydrogen, 
colored according to their metallicity. The gray-shaded 
contours indicate regions of likelihood, based on the pa- 
rameter value ranges listed in column 2 and 3 of Tab. [TJ 
"Typical LAE" models all lie within the dark gray area, 
while "unusual" models lie within the light gray area. 

The remaining models all have rather unrealistic val- 
ues. We have marked with black circles the ones for 
which nevertheless at least one of the following inequal- 
ities is true: nHi, c ] < 1-5 cm -3 , Z < 1.25, oy, c l > 10, or 
P c \ > 0.10; these models are regarded as "extreme, but 
possibly conceivable" , while the remaining are regarded 
as "extreme and probably inconceivable" . 

6. DISCUSSION 

By far, the majority of the models reveal a boost of 
< 1. Furthermore, the ones that do exhibit b > 1 tend 
to have supersolar metallicity, very high cloud Hi den- 
sities (~ 5-10 x typical values), very low velocity fields 
(both V ut and oy iC i are < 10 kms" 1 ), virtually empty 
intercloud media, and a high fraction (> 80%) of photons 
born in the ICM. 

Of the models labeled "extreme, but possibly conceiv- 
able" (indicated in Fig. [15] by the black circles) , twelve 
has resulted in a boost of b > 1.5. Taking a closer look 
at these model reveals that they all lie close to the (ad- 
mittedly somewhat arbitrary) threshold between "con- 
ceivable" and "inconceivable" . Moreover, although out- 
flow velocity was not considered as a defining threshold, 
all except one have V ou t ~ 1 krns" 1 . One model has 
V ut = 17 krns -1 , but has slightly supersolar metallicity 
(Z = I.IZq), low velocity dispersion (ay iC i = 14 krns" 1 ), 
low emission cloud-correlation (P c \ = 0.15), narrow in- 
trinsic emission line (cr om = 15 krns -1 ), high cloud den- 
sities (nni.ci — 1.5 cm" 3 ) while virtually empty ICM 
(n H i,iCM - 10~ 10 cm" 3 ). 

The twelve models were all resimulated with 10 5 pho- 
tons. These rcsimulations reduced the value of the boost 
somewhat in all but one case (where it increased by a few 
percent). Furthermore, inspecting the emerging spectra 
(Fig. I16p shows that all models exhibit very narrow emis- 
sion lines. In order for the Neufeld scenario to work, 
the Lya photons by definition are not allowed to scat- 
ter much, lest they would diffuse too much in frequency, 
eventually being so far from the line center that they 
would be able to penetrate the clouds. The resimulated 
lines have, on average, a FWHM of 110 krns -1 , and the 
broadest line has a width of only 150 krns -1 . This is 
much narrower than observed Lya lines, which typically 
are many hundreds o f krns -1 broad, and rarely below 
250-3 00 krns- 1 (e g fRhoads et all [200l iFvnbo et al.1 
[20031: iHTetall 12001 IShimasaku et al J 12006ft . As a fur- 
ther consideration, in spite of the vanishing outflow ve- 
locities, almost all of the spectra lack the prominent 
double-peaked feature of Lya radiation escaping a static 
medium. 

Additionally, if a given "too high"-EW LAE is sup- 
posed to be caused by a multiphase medium preferen- 
tially absorbing the continuum, rather than some effect 
increasing the Lya flux, then that LAE should also ex- 
hibit a certain reddening of the continuum. In general, 
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Fig. 15. — "Likelihood map" of the obtained boosts for models with realistic and sub-realistic parameter values. The boost is shown, 
arbitrarily, as a function of average neutral hydrogen column density (ATjjj), as measured from the center to the surface of a galax y. Each 
point represents a galaxy, the physical parameters of which are a realization of random values in the ranges discussed in Sees. l4.1l through 
14.101 and summarized in Tab. [T] The metallicity of a galaxy is depicted by the color of the associated point. Models of "typical" and 
"unusual" parameters all lie within the dark and light gray regions, respectively. The rest of the models have rather unrealistic parameters, 
but the ones that have either c j < 1.5 cm" 3 , Z < 1.25, <rv,cl > 10 kms -1 , or P c \ > 0.10 arc marked by black circles. 
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Fig. 16. — Average spectrum (solid line) and lcr region (gray- 
shaded area) of the marginally realistic models giving a boost larger 
than b = 1.5. With an average (maximum) FWHM of 110 (150) 
kms" 1 . these lines are much narrower than typical observed high- 
EW LAEs. 

however, the opposite seems to be the ca se: In a large 
sampl e of COSMOS LBGs and LAEs, iMallerv et all 
(|2012f ) found a clear anti-correlation between the ratio 



of SFRs calculated from the Lya and from SED fitting, 
re spectively, and E(B — V). Similar results were f ound 
bv lShaplev et all ()2003[ ) and lPentericci et al l (120091) . On 
the scale of individual galaxies. lAtek et al.l ()200cl ) found 
in at least four out of six nearby galaxies that the regions 
of high EW correlate with the regions of low E(B — V). 
Most of the high EW regions is diffuse Lya emission, 
which is probably Lya photons produced in the star- 
forming regions being scattered toward the observer. 

6.1. Alternative scenarios 

Based on the results presented above, it seems that 
the Neufeld scenario will allow for an increased Lya EW 
under very special circumstances only, and thus serves 
as an very unlikely explanation for the observed high 
Lya EWs. The reports of high EWs nevertheless exist 
and warrants an explanation. In the following we hence 
discuss miscellaneous alternatives. 

6.1.1. Measuring errors 

Firstly, it must be remembered that the mere re- 
port of an enhanced EW not necessarily implies an en- 
hanced EW. For faint continua, determining the EW 
is associated with considerable error, and most of pro- 
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claimed high-EWs have huge error bars. For instance, 
iHenrv et all (|2010j) found six sources with W cul > 240 
A, but caution that the uncertainties in such high EWs 
typically exceed 100 A. 

6.1.2. Stellar population 

The the oretical upper l imit for the Lya EW of 240 A 
assumes a lSalpete"n (]1955[ ) IMF at Solar metallicity. For 
decreasing metallicites, the main sequence of a stellar 
population is shifted to the blue, resulting in a slower 
decline of the io nizing photon production, and thus a 
higher Lya EW jSch aererl [20031) . Nevertheless, even for 
Z ~ 0.02Z Q , the EW of a starburst is "only" 350 A, i.e. 
corresponds to an interpreted boost of ~ 1.5. 

For even lower values of Z very high values of W[ n t 
may be reached. An extreme case of a low-mctallicity 
population would be one consisting of PopIII s tars (e.g. 
iTumlinson et all 120031: iZackrisson et all 120111) . While 
still not observationally confirmed, stars born from (close 
to) metal-free gas are expected to be able to reach very 
high masses. Being particularly luminous in the Lyman 
continuum, such a population of stars will be able to emit 
a remarkable fraction of their photons in Lya. Whereas 
a "normal" population emits 6-7% of its bolometric lu- 
minosity in Lya, 10-40% can be reached for sufficiently 
low metallicities, with a resulting W- m t of several thou- 
sands of Angstrom ([Bromm et al J 1200 U: iSchaererl [20021: 
iRaiter et alj l20l"oT ). However, at metallicities this low, 
no dust should yet have formed, and thus no reddening 
should be observed. Moreover, with such a hard spec- 
trum, the EW of the Hen Ha line at 1640 A is expected to 
be noticable. In some of the high Lya-EW cases, this was 
speci fically looked for, but not detected ([Dawson et al.l 
l200l . 

Shifting one or both of the IMF mass limits toward 
higher masses, or applying a shallower slope, implies a 
relatively larger fraction of high- mass stars, with a harder 
ionizing UV spectrum and thus a larger Wi n t. For in- 
stance, increasing the upper mass limit from 100 to 500 
Mq (corresponding to a rise in average mass from 3.1 to 
3.4 Mq only) for a Z = 5 x 10~ 4 Z© population, strength- 
ens Wi n t from ~ 400 A to ~ 700 A, while further increas- 
ing the lower mass limit from 1 to 500 Mq (corresponding 
to a rise in (M) to 112 Mq) strengthens W int to ~ 900 
A (see Fig. 7 of lSchaererll200l . 

In all cases, the above values assume an instantaneous 
burst of stars; with the early O and B star dying fast, 
the EW of these extreme populations quickly decline, 
reaching more modest values after only a few Myr, and 
would thus have to be observed in a special period of 
their lives. For continuous star formation, high EWs can 
be found at any point in time, although the values will 
always be lower than the initial EWs of the instantaneous 
bursts. For Z = 5x 10~ 4 Zq, Wi nt increases from ~ 180 
A to - 240 (~ 500) A for a mass range 1-500 (50-500) 
Mq. This can be compared to W ln t ~ 100 A or less for 
Z > 0.02 Zq. 

6.1.3. Delayed escape of Lya 

In principle, for a short-lived starburst the Lya pho- 
tons — having their path length out of the galaxy in- 
creased — could be observed after the continuum has 



already faded, resulting in a higher EW. This was in- 
vestigated for a homo geneous sp here of dust-free gas by 
IRov et al.l (|2010l ) and Ku et alT(|2011l ), who found that, 
depending on the column densities, the Lya radiation 
may be appreciably delayed. When dumpiness is intro- 
duced, we find that the escape time decreases signifi- 
cantly. For instance, if all the gas in the fiducial model is 
distributed homogeneously, the Lya photons escape with 
a typical delay of t csc ^100 kyr if there were no dust. 
In contrast, with the clumpy structure, t csc ~ 20 kyr 
(not considering the photons that escape directly with- 
out interacting with any clouds). When dust is added, 
t csc is reduced further, since the photons that have the 
longest escape time are the ones that are most prone to 
absorption. 

For larger r ga i, t csc increases. However, even among 
the models with r ga i = 10 kpc, few have t csc > 100 kyr. 
Compared to the typical lifetime of the O and B stars of 
several Myr, a stellar continuum should still be visible as 
the Lya radiation reaches its maximum. 

Solving simultanenously the RT for Lya and ionizing 
UV radiaton in a sp herically symmetric model galaxy, 
([Yaiima et al.ll2012d ) found that most of the Lya pho- 
tons could remain trapped until the galaxy has been fully 
ionized, at which point the "old" and the "new" Lya pho- 
tons are released together, resulting in Lya EWs of up 
to ~ 1000 A. It is unclear, however, how applicable these 
results are to more realistically simulated galaxies, where 
low-density regions are ionized more rapidly, allowing the 
Lya photons to leak in those directions. 

In a cosmol ogically simulated galaxy of appro ximate 
radius 25 kpc, lLaursen fc Sommer-Larsenl (j2007| ) found 
a typical path length of the Lya photons of ~ 40 kpc, i.e. 
a delay with respect to the continuum of ~ 50 kyr. This 
simulation, however, did not include ionizing UV RT. 

6.1.4. AGN activity 

If a galaxy contains an AGN, the Lya EW may be 
much higher than the 240 A. The line widths of these 
object are so broad, however, that such ob jects can usu- 
ally b e ex cluded, as in the c ase of e.g. LRhoads et al.1 
(|200l and lDawson eTaTl ([200l . Even type II AGNs, i.e. 
objects where the broad-line region is not visible, typi- 
cally have larger line widths than the observed high-EW 
objects. Furthermore, if the majority of the high-EW 
systems were due to AGN activity, i t would be i ncon- 
sistent with the X- ray emission (e.g. iWang eFaLl [20P: 
iGawiser et al]|200g) . 

6.1.5. Viewing angle 

All the models studied in this work are spherically sym- 
metric. For a flattened system, photons escape more eas- 
ily face-on than edge-on. But since Lya photons scatter 
and thus wander around in the galaxy, their probabil- 
ity of escaping from the face before reaching the edge is 
increased, and thus a higher EW will be observed face- 
on. Without the information of the system's morphology, 
this can be interpreted as an EW boost. That this is true 
even for a homogeneous, dust-free ISM is seen in Fig. 1171 
where the impact on the EW of gradually flattening a 
sphere of total Hi column density N-hj = 2 x 10 20 cm -2 
(corresponding to a DLA) is shown. 

This effect is of course observationally impractica- 
ble, except statistically, but has also been investigated 
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Fig. 17. — EW "boost" b when observing the system face-on 
(solid) and edge-on (dotted) as a function of flattening g = 1 — c/a, 
where a and c are the semi-major and -minor axes of the sys- 
tem (shown schematically as small ellipses in the top of the plot). 
Results are shown for both dust-free (red) and dusty (green, 
Z = O.IZq) gas. In both cases the gas is homogeneous, and the 
column density is ./Vjjj = 2 X 10 20 cm -2 , as measured along the 
c-axis. For reference, the M ilky Way has a flatness of g ~ 0.95-0.99 
UBinnev fc M crrificld 1998). 



in mor e realistic galaxy simulations. In lLaursen et al.l 
(|2009bD . Lya RT was carried out in a sample of nine 
randomly oriented galaxies extracted from a fully cos- 
mological simulation and resimulated at high resolution. 
The ratio of the fluxes escaping in the most luminous to 
the least luminous direction, roughly corresponding to 
what would be interpreted as an EW boost, was in the 
range 1.5-4 if integrating over the whole galaxy, and ~ 10 
if looking at the region of maximum surface brightness 
only. From Fig. |T7l this corresponds roughly to "boosts" 
of 1.1—1.5 and ~ 2, res pectively. Similar values was found 
by lZheng et al.l ()2010[ ) (a ra tio of seven for the ce ntral re- 
gion of a single galaxy) and iBarnes et al.l (j2011l ) (a ratio 
of 1.7-3 for va rious implementation s of galactic winds). 
More recently, lYaiima et al.l (|2012bD found no preferred 
direction of escape for a single, cosmologically simulated 
galaxy at z = 3.1, but a strong dependence on viewing 
angle at z = when the gal axy had settled into a d isk, 
similar to what was found by iVerhamme et al.l ()2012l ) for 
an isolated disk galaxy. 

6.1.6. Cooling radiation 

At high redshifts, a large fraction of galaxies may be 
expected to be in the process of forming, i.e. to accrete 
gas from the surrounding medium, resulting in collision- 
ally excited Hi leading to Lya emission. Since the im- 
plied densities are rather high, the timescale for cooling 
is generally smaller than the dynamical timescale, and 
thus t he infalling ga s sho uld be relatively cold (T ~ 10 4 
K, e.g. lFall fc ReelH985l: iHaiman fe Reesll200lD . Conse- 
quently, ro ughly 50% of the e nergy is emitted in the Lya 
line alone (jFardal et al.ll2001l) . 

For relatively massive galaxies, iGoerdt et al.1 (j2010l ) 
found at z = 2.5 that cold gas accreting can give 
rise to significant Lya emission of up to a few times 
10~ 17 erg s _1 cm -2 arcscc -2 on scales of 10-100 kpc, 
which is comparable to typical line strengths. Such 
gravitational cooling not only gives rise to Lya radia- 
tion but also to a continuum, although relatively it is 



much fainter than the stellar continuum. iDiikstral ([2009D 
calculated analytically the EW of the cooling radiation 
from a collapsing halo and found that, depending on 
where the continuum is measured, the intrinsic EW of 
the cooling radiation is > 950 A. Numerically, a few 
studies have addressed this issue and find that in gen- 
eral the fraction of the Lya radiation that stems from 
gravitati onal cooling increases with z, fr om ~ 10% at 
z = 3.6 (lLaursen fc Sommer-Larsen1 l2007D. to ~ 16-18% 
at z = 5.7^6.5 (iDaval et al.||2010l), to ~ 50-60% a t z = 7- 
8 (jJensen et al.ll2012D . Similarly. lYaiima et ail (|2012al) 
found that cooling radiation becomes comparable to re- 
combination Lya for z > 6. Obviously, large deviations 
from these values exist, both numerically and in real- 
ity, as in the case of Lya blobs where in some cases 
no stellar continuum is seen at all; a liable explanation 
for these objects is co ld accretion (|Nilsson et al.l 120061 : 
iDiikstra fc Loebll2009l ) 

However, while the photons of stellar origin, both Lya 
and continuum, tend to be born in regions of compara- 
tively high densities of both dust and neutral gas, thus 
being more prone to dust absorption, the cooling pho- 
tons may escape more or less freely. Thus, cooling radia- 
tion could dominate even at quite low redshifts, although 
the actual fraction is difficult to assess observationally as 
long as the escape fraction of stellarly produced photons 
is still associated with so large uncertainties. 

6.1.7. IGM extinction 

At high redshifts, neutral hydrogen in the IGM scat- 
ter a significant fraction of the photons blueward of the 
Lya line out of the line of sight. This affects both 
the continuum and the line itself, and is sometimes ac- 
counted for using a model of the IG M absorption. For 
instance, iMalhotra fc Rhoadsl (|2002t) use the prescription 
of IMadaul (11991 " and found that the corrections for the 
broadband and the narrowband approximately cancel, 
leaving the EW unaffected. This assumes that the Lya 
line is symmetric about the line center; if the red peak is 
enhanced relative to the blue peak, as would be the case 
for an expanding medium, the IGM transmits more Lya, 
resulting in the erroneous interpretation that the EW is 
larger than in reality. Moreover, as the galaxies detected 
in Lya at high redshifts may be biased toward being the 
ones that have made themselves visible by ionizing a bub- 
ble around them in an otherwise partly neutral IGM, the 
line may be less affected by the IGM than the continuum, 
which again would be interpreted as an EW boost. 

Even if the correlation of the state of the IGM with 
the sources is taken into account, large fluctuations 
in the neutral fraction of hydrogen exis ts, which may 
also be a possible source of e rror (e.g. lLaursen et al.l 
120111 : iJeeson-Daniel et al.l l2012h . Using a simple model 
for LAEs at high redshifts bu t incor porating dumpi- 
ness of the IGM, lHaves et al.l (|2006l ) calculate distri- 
bution functions for observed EWs, for various narrow- 
band/broadband techniques. They found that, depend- 
ing on the redshift, the median W^bs may easily lie ap- 
preciatively above W em , and that for instance at z = 5.7 
with narrowband and broadband both centered on the 
Lya line, 20% of the EWs will be boosted by a factor of 
2 or higher. 

6.1.8. Star formation stochasticity 
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When calculating the resulting EW from a dis- 
tribution of stellar masses, an infinite number of 
stars are u sually assumed. Using the numerica l 
code SLUG (iFumagalli et alJl2?niH d a Silva et alJ lifflllh. 
iForero-Romero fc Diikstral (|2012f ) found that for low- 
SFRs galaxies, stochasticity alone results in a distribu- 
tion of EWs from Wo/4 to 3 Wo, where Wo is the mean 
EW. However, since this effect is minimal for SFRs > 1 
M Q yr _1 , which are hardly visible at high redshifts, it 
cannot serve as an explanation for most of the observed 
high-EW galaxies. 

6.1.9. Inhomogeneous escape 

For resolved galaxies, regions with no star formation 
may scatter Lya from distant regions toward the ob- 
server, such_that_^ese_particular regions exhibit extreme 
EWs (|Haves et al J 12007( 1 . Integrated over the source, 
however, Lya-overluminous regions cancel out with un- 
derluminous regions (modulo other effects affecting the 
RT). 

7. SUMMARY AND CONCLUSION 

We have carried out a comprehensive set of Lya and 
FUV continuum radiative transfer calculations, system- 
atically varying the significant physical parameters, first 
separately and then in unison. The aim was to inves- 
tigate whether a multiphase interstellar medium is ca- 
pable of preferentially absorbing the continuum so as to 
enhance the observed Lya equivalent width. The moti- 
vation for this study are the numerous observations of 
Lya emitting galaxies exhibiting EWs larger than what 
is theoretically expected, together with the previously 
proposed explanation of s uch a "boost" being a conse- 
quen ce of a clumpy ISM (INeufeldl fl99lt E anscn & Oh 
:. In the suggested model, a galaxy consists of a 
number of clouds, physically corresponding to the ISM 
phase termed the "warm neutral medium" , dispersed in 
a low-density medium corresponding to the phase called 
the "hot ionized medium" . 

We find that, while indeed it is possible to construct a 
model galaxy with the said ability, the physical proper- 
ties needed are extremely unlikely to be found in a real 
galaxy. In order for a galaxy to enhance its intrinsic Lya 
EW, a large number of criteria must be met: 

1. The metallicity must be high (> Solar), 

2. the density of the WNM must be very high (> 5 
times typical values), 

3. the density of the HIM and its ionization fraction 
must be very low, so that the total density of neu- 
tral hydrogen in this phase does not exceed ~ 10~ 7 
cm~ 3 , 

4. the bulk (> 80-90%) of the Lya and FUV photons 
must originate from regions where the stars have 



blown away the neutral gas from which they were 
born, 

5. the galaxy must have virtually no outflows (K» u t < 
10 kms' 1 ), and 

6. the velocity dispersion of the gas clouds must be 
very low (< 10 kms -1 ). 

In particular the latter two points are important; the 
whole point of the proposed multiphase model is that 
the paths of the Lya photons are confined to a dust- and 
gasless medium, never penetrating the dusty clouds but 
instead scattering off of their surfaces. This mechanism 
works because of the resonance nature of Lya scattering. 
However, as soon as the clouds have a small velocity, es- 
pecially random motions, the Lya photons are shifted 
out of resonance in the reference frame of the clouds, al- 
lowing them to penetrate much farther into the clouds 
and thus exposing them to a higher column density of 
dust. Note that in principle it is possible to have an en- 
hanced EW for larger velocity fields, but this will require 
even more extreme conditions regarding the other param- 
eters. With regards to the required densities of neutral 
hydrogen, these are more characteristic of the so-called 
cold neutral medium or molecular clouds. Since the vol- 
ume occupied by t hese phases are much smaller than that 
of the WNM (e.g. iFerrierd 12001 . the covering factor of 
these phases is small. Moreover, they will typically be 
embedded inside the WNM, affecting only those photons 
that were already trapped in the clouds. In addition to 
the points mentioned above, it must also be remembered 
that the intrinsic Lya EW of a galaxy exhibiting a high 
observed EW must be high, i.e. the stellar population 
must be very young (~a few Myr); otherwise an even 
higher boost would be needed to bring the EW above 
the theoretical limit. 

In conclusion, we consider the Ncufcld model to be an 
extremely unlikely reason for the observed high EWs. 
We have discussed a number of other possible explana- 
tion (Sec. 16. 1| . of which the most probable, in our opin- 
ion, is excess of Lya radiation from cold accretion (may 
contribute ~ 0.1-1 times the stellar Lya, depending 
on redshift, but may escape much more easily), and/or 
anisotropic escape of Lya (the Lya flux may vary in dif- 
ferent directions by a factor ~ 1 10, depending on irreg- 
ularity of the galaxy and aperture of the detector. 

We are grateful to Jens Knude, Sangeeta Malhotra, 
and Andrei Mesinger for valuable discussions about the 
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fundings from the Villum Foundation. The simulations 
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APPENDIX 

EXTENDING AND TESTING MoCaLaTA 

Implementing continuum radiation 

The original version of MoCaLaTA included only the RT of Lya, i.e. the initial frequency probability distribution 
of photons followed a Voigt profile with the Gaussian core given by the temperature/turbulence of the gas. To calculate 
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EWs, the wavelength range AA transferred though the medium must be extended so as to include the continuum. 
For a given EW this is readily accomplished by first determining whether an emitted photon is a line rather than a 
continuum photon, with a probability P(line) cx W; n t/AA, and then determining the exact wavelength of the photon 
depending on the intrinsic line profile or the slope of the continuum (in this study taken to be flat in v). For a broad 
wavelength region around Lya, numerically there is no difference in the RT of the two kinds of photons. 

Building the grid 

In order to construct spherical clouds, we make use of adaptive mesh refinement (AMR), where cells may be refined 
into eight subcells. As the clouds themselves are uniform, high resolution is only needed at their surfaces, not their 
interior. To build the AMR grid, we use the following approach: 

Given a number N c \ of clouds and a distribution of sizes specified by a minimum and maximum radius r min and 
r mM , as well as a size distribution power law index j3, cloud centers are placed randomly in a sphere of radius r ga i, 
such that no clouds overlap. The clouds' surfaces are then marked by a number of auxiliary particles, typically ~ 10 3 . 
Subsequently, the grid is constructed starting from a mother grid of base resolution 2 3 , recursively refining any cell 
that contains more than one particle. Finally, each cell is labeled either a "cloud" cell or an "ICM" cell, depending 
on the distance to the nearest cloud center, and the physical parameters then assigned accordingly. Note that all cells 
belonging to a given cloud share the same bulk velocity, such that there are no velocity gradients inside a cloud. 

The auxiliary particles can be positioned in different ways so as to constitute the cloud surface; the most straight- 
forward way is to place them at a random position at a distance r c i from the center of the cloud. However, it turns 
out that in general some particles end up so close to each other, that in order to refine the cells sufficiently, hundreds 
of refinement levels may needed. Distributing the particles evenly on the surface ensures that rough ly the same level 
of refi nement is needed everywhere. Various algorithms exist for this geometrical exercise (see, e.g.. iSaff k, Kuiilaarsl 
[19971) . 

Altering the acceleration scheme 

To speed up calculations, MoCaLaTA makes use of various acceleration schemes. The most efficient is the core- 
skipping scheme: Since scatterings of photons near the line center are usually associated with negligible spatial move- 
ment, while it typically takes of the order 10 5 scatterings to "escape" the core, these scatterings can be skipped by 
drawing the random velocity of the scattering atom from a centrally-truncated Gaussian, favoring fast moving atoms 
that can give the photon a large Dopplcr shift. In MoCaLaTA, the critical wavelength defining when the acceleration 
scheme can be invoked is dependent on the product aro, such that for higher Hi column densities and lower temper- 
atures, a larger critical value can be used, allowing for a larger speed-up. The value of aro is specific to a given cell, 
and is calculated from the center of the cell to the face. 

Since the code has previously been used for galaxies extracted from cosmological simulations, cells are usually 
surrounded by other cells of similar physical conditions, i.e. there are no large gradients. In the present, idealized 
simulations, however, where densities may go from zero to some extreme value over the course of a single cell step, 
this poses a problem for photons close to the line center scattering on the surface of a cloud. 

In reality, the photon typically escapes after a handful of scatterings, changing its frequency only of the order 1.5 
Doppler widths. With the acceleration scheme, the first scattering pushes the photon artificially far out in the wing, 
rendering the optical depth of hydrogen significantly reduced, such that the photon is able to penetrate the cloud (or 
the next cloud) too much, with a corresponding higher probability of being absorbed by dust. 

We solve this issue by calculating otq not from the cell center, but from the current position of the photon to the 
face of the cell. Unfortunately, the consequence is a quite less efficient acceleration. 

Testing the code 

In the original Ncufeld scenario, that is with no velocity fields etc., the physical parameters most important to the 
RT are Nq and e c , i.e. the average number of clouds with which a Lya photon emitted from the center of a spherical 
conglomeration of clouds interacts, and the probability per cloud interac tion that a photon (Lya or continuum) 
is absorbed. No analytical solution exist for these quantities, but IHO06I found that they could be well fitted by 
Ao = / c 2 + f/c (in the absence of absorption) and e c ~ 3e^ 9 / (1 + 2eJ 2 ), where is the probability of absorption per 
[gas or dust] interaction for an incident photon lying at a;, Dopplcr widths from the line center. The results are shown 
in Fig. QB 
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